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Dome-shaped magnetic phase diagram of thermoelectric layered cobaltites 
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Using muon spin spectroscopy we have found that, for both Na^CoCh (0.6 < x < 0.9) and 3- and 
4- layer cobaltites, a common low temperature magnetic state (which in some cases is manifest as an 
incommensurate spin density wave) forms in the C0O2 planes. Here we summarize those results and 
report a dome-shaped relation between the transition temperature into the low-T magnetic state 
and the composition x for Na^CoC^ and/or the high-temperature asymptotic limit of thermopower 
in the more complex 3- and 4-layer cobaltites. This behavior is explained using the Hubbard model 
on two-dimensional triangular lattice in the C0O2 plane. 

PACS numbers: 76.75.-K, 75.30.Fv, 72.15.Jf, 75.30.Kz 



Although the widespread current interest in the lay- 
ered cobaltites 0, Q was originally due mainly to their 
unique combination of high thermopower S with metal- 
lic transport properties, 0, 0, S IM S E3, EH E3l 
which makes them one of the most promising systems for 
power applications, we have shown that they also display 
interesting and complex magnetic orderings, directly cor- 
related with the enhanced thermopower. The richness of 
behavior of the layered cobaltites is due to their intrin- 
sic structure, namely: electrically active triangular plane 
layers of C0O2, which are separated by a variety of in- 
termediate structures; their relatively strong electronic 
correlations; and the fact that the structures between 
the C0O2 layers can be modified in a variety of ways to 
vary their dimensionality, ionic states, carrier doping in 
the C0O2 planes and the relevant interaction strengths. 

In order to elucidate the magnetism in the C0O2 
planes and the mechanism of the good thermoelectric 
properties, we have carried out positive muon spin ro- 
tation and relaxation (/i + SR) experiments on the lay- 
ered cobaltites. As a result, we found the transi- 
tion from a high-temperature paramagnetic to a low- 
temperature commensurate or incommensurate spin den- 
sity wave (C- or IC-SDW) state for [Ca 2 Co0 3 ]^|,ICo0 2 ] 
below -TOO K,[H 14] Na . 75 CoO2 at 22 KJ13 and 
[CaaCo^gCua/sO^JafCoOa] below ~200 K,[l|] where 
RS denotes the rocksalt-type subsystem. The common 
magnetic ordering is thus not always that of a classic 
frustrated AF (triangular planar arrangement of spins) 
system but rather a variety of states including disordered 
AF and (most common) an incommensurate spin density 
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wave (IC-SDW), n arising in parallel with other effects in 
the triangular lattice such as effective mass enhancement 
and the enhanced S in the Na^CoC^ case. The more 
complex cobaltites exhibit, in addition, other magnetic 
orderings and even higher S. 

Nevertheless, the transition temperature was scat- 
tered in the wide temperature range (22-200 K), and 
there was no clear relationship between Tsdw an d struc- 
tural properties and/or carrier concentration of the lay- 
ered cobaltites. We back therefore to the basic system, 
Na x Co02, to investigate the dependence of Tsdw on x 
by /i + SR. Also, for the multilayer systems, we clarify 
the magnitude of Tsdw as a function of the Co valence 
varied by the change in composition of the rocksalt-type 
subsystem. [H, 0, UJJ For the Na^OoC^ case, we show 
that a model based on the Hubbard model for a trian- 
gular lattice of S=l/2 Co ions with occupancy given by 
x (correlated with band filling) [ItI ITsj is able to explain 
the dome shaped magnetic phase diagram shown below. 
Similar behavior is observed for the multilayer systems, 
where the increased two dimensionality due to the extra 
layers leads to higher Tsdw for the IC-SDW state and is 
correlated with higher S. 

Single-crystal platelets of Na^CoC^ were prepared by a 
flux method. [19j Polycrystalline 3- and 4-layer cobaltites 
listed in Table 1 were synthesized by a solid state re- 
action technique H, 0, ^| or a reactive templated 
grain growth technique. [2(| The /z + SR experiments were 
performed on the M15 and M20 surface muon beam 
lines at TRIUMF. The experimental setup is described 
in elsewhere. Illll^ 

Figure 1 shows zero-field (ZF-)/x + SR time spectra at 
1.8 K for single crystal platelets of Nao.gCoC>2. The top 
spectrum was obtained with the initial /i + spin direction 
5^(0) perpendicular to the c-axis and the bottom one 
with Sfj, (0) parallel to c. A clear oscillation due to quasi- 
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FIG. 1: ZF-/i + SR time spectra of single crystal platelets of 
Nao.gCo02 at 1.8 K. The configurations of the sample and 
the initial muon spin direction (0) are (top) (0) _L c and 
(bottom) S M (0) |j c. The bottom spectrum is offset below to 
be seen clearly. 



static internal fields is observed only for S^(0) perpen- 
dicular to c. The muon signal is fitted best by a zeroth- 
order Bessel function of the first kind, Jo(w^t), which 
describes the evolution of the muon polarization in an 
IC-SDW field distribution; |2l|.[22| we therefore conclude 
that Nao.gCo02 undergoes a magnetic transition from a 
paramagnetic state to an IC-SDW state at Tsdw=19 K. 
The absence of a clear oscillation for the S^(0) || c case 
indicates that the internal field Hint is roughly parallel 
to the c-axis. Because of the strong anisotropy, the IC- 
SDW is most likely to propagate in the c plane (i.e. in 
the C0O2 plane), with oscillating moments directed along 
the c-axis. 

Figure 2 shows the temperature dependence of the 
muon precession frequency Vn =lo^/2tt for the single crys- 
tal platelets of Nao.gCoC^. Here is the order param- 
eter of the transition and its T dependence is well de- 
scribed by the BCS weak coupling expression for such 
order parameters, as expected for the IC-SDW state. |23l] 
The magnetic phase diagram (Fig. 3) of Na^OoC^ can 
thus be sketched from the /x + SR results for polycrys- 
talline samples with x=0.65 and 0.75 [15J and the present 
x=0.6 and 0.9 crystals. Recent compositional and chemi- 
cal titration analyses indicated that the oxygen deficiency 
5 in Na- £ Co02- ( 5 is negligibly small even for the £=0.9 
sample. [T3 Hence, the average Co valence can be directly 
calculated from x. As x increases from 0.6 (i.e. the Co 
valence decreases from 3.4), the magnitude of Tsdw in- 
creases up to around x=0.8, then decreases with further 
increasing x. As a result, we obtain the dome-shaped 
relationship of Figure 3, between Tsdw and x, viz. the 
Co valence for Na^CoC^. 

The x=l end member of Na^CoC^, that is, a fully 



FIG. 2: Temperature dependence of the muon precession 
frequency in Nao.gCo02- The solid line represents the 
temperature dependence of the BCS gap energy. 
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FIG. 3: Phase diagram of Naj,Co02 determined by the /x + SR 
experiments. Solid and open circles represent the present re- 
sults on single crystals and polycrystalline samples, respec- 
tively. The point at x=l is extrapolated from the data on the 
related compound LiCoQ9. |24| 



occupied NaCoC>2 phase, cannot be prepared by con- 
ventional solid state reaction and/or flux techniques, 
whereas the related compound LIC0O2 is easily obtained. 
The structure of LiCoC>2 is isomorphous with a-NaFcC>2 
(i?3m), almost the same as that of Na^CoOa, and it is re- 
ported to be diamagnetic down to 4.2 K.j2J] This behav- 
ior is expected, because the Co 3+ ions in LiCoC>2 are in a 
low-spin state t\ g , as for Na x Co02- Therefore, Na^OoC^ 
is also expected to lack magnetically ordered states. 
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The occupancy of Co 4+ spins (S=l/2) in the two- 
dimensional triangular lattice (2DTL) increases with de- 
creasing x. Thus the other end member, NaoCo02, would 
be a half filled 2DTL. In other words, every lattice site is 
occupied by an S=l/2 spin. The Hubbard model within 
a mean field approximation can be used for explaining 
the magnetism of such a system, with the Hamiltonian 
H given by [H|l3 



H 



■ C itr C j<? 
<ij>cr 



U J2 



(1) 



where c\ a {cj a ) creates (destroys) an electron with spin a 
is the number operator, t is the 



on site i, ru 



- J 



nearest-neighbor hopping amplitude and U is the Hub- 
bard on-site repulsion. The electron filling n is defined as 
n=(l/2N)'^2 i rii, where N is the total number of sites. 

At T=0 and n=0.5 (i.e., NaoCo02), as U increases 
from 0, the system is a paramagnetic metal up to 
U/t=3. 97 due to geometrical frustration, then changes 
into a metal with a spiral IC-SDW, and then at U/t=5.27 
a first-order metal-insulator transition occurs. Il7j The 
lack of magnetic transitions for Na x Co02 with x=0.6 
and 0.65 suggests that U/t <3.97. This means that 
Na x Co02 is unlikely to be a typical strongly correlated 
electron system, because U ^> t for such a system. The 
calculations^^ also predict that, as n increases from 
0, the magnitude of U jt at the boundary between the 
paramagnetic and SDW phases decreases, with increas- 
ing slope (d(U/t)/dn) upton=0.75. Even for U/t=0, the 
SDW phase is stable at n=0.75. The value of U /t then in- 
creases with further increasing n, with decreasing slope. 
Therefore, the dome-shaped phase diagram of Fig. 3 is 
qualitatively explained by the calculations, although the 
measured maximum of the dome is located around cc=0.8 
(i.e. n=0.9). This is likely due to the simple band 
structure assumed in the above calculation, while cal- 
culations for Na^CoC^ suggest more complicated band 
structure. |25| 

For the related compounds [Ca2CoC>3]^ ^[CoC^] and 
[Ca2Co4/ 3 Cu2/304]^| 2 [Co02], (i.e. 3-layer and 4-layer 
cobaltites, where RS denotes the rocksalt-type susbsys- 
tem), the IC-SDW transition was also observed as a com- 
mon behavior in the C0O2 planes. [T3I IbH Il6| We have 
studied by /i + SR the dependence of the IC-SDW transi- 
tion on the Co valence for the variety of layered cobaltites 
shown in Table 1. It should be noted that the magnetic 
susceptibility x(^) curves for the 3- and 4-layer samples 
lacked a marked change at either the or Tsdw de- 

tected by ^ + SR. 

There are two Co sites in these layered cobaltite lat- 
tices, one in the rocksalt-type subsystem and the other 
in the C0O2 plane. Therefore it is difficult to determine 
the Co valence in the C0O2 plane by x measurements or 
chemical titration techniques. Since the Hall coefficients 
i?H of these cobaltites are reported to depend strongly on 
temperature, 0,^1 it is a ls° difficult to estimate the car- 
rier concentration in the C0O2 plane. On the other hand, 
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FIG. 4: The relationship between Tsdw and thermopower 
S at 300 K. Solid circles represent the data for the 
cobaltites with a triple rocksalt-type subsystem, open di- 
amonds a quadruple rocksalt-type subsystem and crossed 
squares Na^OoO^. 



the magnitude of 5 in the layered cobaltites is almost in- 
dependent of T above 200 K. 0, S H EIH El There- 
fore, we use the value of S at 300 K as an indicator of the 
Co valence in the C0O2 plane for these cobaltites. The 
modified Heikes formula gives 5 at the high-temperature 
limit (St^oc) as ^ 



Si 



IS 



-^ln(^ 

e \04 1 - V 



(2) 



where is the Boltzmann constant, e is the elementary 
charge and g 3 and g± are the numbers of the spin con- 
figurations of Co 3+ and Co 4+ ions, respectively, and y 
is the ratio Co 4 +/(Co 3+ +Co 4+ ). Since both Co 3 + and 
Co 4+ are in the low-spin state (t® g and t\ g ) and 53=1 
and <?4=6, we can convert 5(300 K) to y, assuming that 
5(300 K)=5 T ^oo. 

Figure 4 shows T|d W as a function of 5(300 K) for 
the all cobaltites listed in Table 1. A clear dome-shaped 
relation is thus observed for all the available cobaltites 
with a variable number of layers (N) between the two 
adjacent C0O2 planes. As N increases from 1, the 
Tgp W — vs. — 5(300 K) curve shifts towards higher tem- 
perature. This is due to the increased two-dimensionality 
induced by the increase in the interlayer distance be- 
tween C0O2 planes. Also, the large observed transi- 
tion widths (50-130 K) are consistent with enhanced 
two-dimensionality and resulting spin fluctuations. Phe- 
nomenologically, the phase diagram is very similar to the 
well-known relationship between the superconducting T c 
and the Cu valence in the high-T c cuprates. Actually, 
both SDW and superconducting transitions are induced 
by an intrinsic instability of an electron system; in other 
words, as T decreases, an energy gap appears at Tsdw 
and/or T c to minimize the internal energy for both cases. 
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TABLE I: Parameters of the magnetic transition and thermopower S of several cobaltites; the number of the layers between 
the two adjacent [C0O2] planes (N), the onset and endpoint temperature of the magnetic transition detected by the /i + SR 
experiments and S at 300 K. nd>3.3 means not detected down to 3.3 K. * (**) indicates the sample showing a clear muon 
precession due to the IC-SDW (C-SDW) field in a ZF-/1+SR spectrum. 



cobaltite 


iV 


moil (Tr\ 

TiSw (K) 


Isdw (K) 


S(300 K) (fiVK ) 


[Ca 2 Co0 3 ]o:62 [C0O2] 


3 


100±5 


30±5* 


130 


[Cai. 8 Bi .2CoO 3 ]^ b [CoO2] 


3 


120±20 


70±10* 


140 


[Cai. 8 Y .2CoO 3 ]" [C0O2] 


3 


120±20 


70±10 


140 


[Cai.sSro^CoOaJ^ [C0O2J 


3 


100±10 


45±10* 


118 


[Pb .4Coo.6Ca 2 03]a62 [C0O2] 


3 


130±20 


40±10* 


165 


[Pbo.7Coo.4SriCao.90 3 ]o:58[Co02] 


3 


90±10 


10±10 


130 


[Pb().7Coo. 3 Sr20 3 ]o.56[Co0 2 J 


3 


80±10 


nd>3.3 


120 


[Tlo.76Co . 33 Sr 2 3 ]5! 6 [Co02] 


3 


10±40 


nd>5.1 


90 


[Tli.oiCoo.ii Sri. 88 3 ]£i 6 [Co0 2 ] 


3 


25±5 


nd>1.7 


90 


[Tli .is Coo.oa Sn.reOslSir [C0O2 ] 


3 


nd>2.5 


nd>2.5 


90 


tHgo.4Coo.6Sr 2 O s ]Sl6[C!o0 2 ] 


3 


40±25 


5±5 


100 


[Tio.4Coo.6Ca20 3 ]5| 2 [Co02] 


3 


150±10 


15±5 


150 


[Ca 2 C04/3 C3u a /30 4 ]aS2 [Co0 2 ] 


4 


190±10 


145±10* 


150 


[Ca 2 Bii. 7 Co . 3 04]o l 6o[Co0 2 ] 


4 


60±20 


nd>5 


140 


[Sr 2 Bi 2 4 ]P[Co0 2 ] 


4 


nd>2.3 


nd>2.3 


134 


[Ba 2 Bi 2 04]o:![Co02] 


4 


nd>1.8 


nd>1.8 


100 


Nao.gCo0 2 


1 


19±0.5 


19±3* 


150 


Na .75CoO 2 


1 


22±0.5 


22±r* 


120 


Na .65CoO 2 


1 


nd>2.5 


nd>2.5 


100 


Nao.eCo0 2 


1 


nd>4.3 


nd>4.3 


100 



Therefore, it is reasonable to expect a similar relation- 
ship between transition temperature and carrier concen- 
tration for both the magnetic cobaltites and the super- 
conducting cuprates. 

Furthermore, the average Co valence for the maximum 
Tsdw indicates the optimal filling to induce an SDW tran- 
sition at high temperatures and enhance the effective 
mass of charge carriers through the AF interaction be- 
tween spins. In other words, this dome relation provides 
important guidance in the search for improved thermo- 
electric properties of the layered cobaltites. 

We thank S.R. Kreitzman, B. Hitti, D.J. Arseneau, 
Y. Seno, K. Dohmae, C. Xia, H. Nozaki, H. Hazama, J. 



Chakhalian, D. Liu, A. I.-Najafabadi and S. D. LaRoy 
for help with the /i + SR experiments. Also, we appreci- 
ate R. Asahi, U. Mizutani, H. Ikuta, T. Takcuchi and 
K. Machida for discussions. This work was supported 
at Toyota CRDL by joint research and development 
with International Center for Environmental Technology 
Transfer in 2002-2004, commissioned by the Ministry of 
Economy Trade and Industry of Japan, at UBC by the 
Canadian Institute for Advanced Research, the Natural 
Sciences and Engineering Research Council of Canada, 
and at TRIUMF by the National Research Council of 
Canada. 



[1] K. Takada et al. Nature 422, 53 (2003). 

[2] R. E. Schaak et al., Nature 424, 527 (2003). 

[3] J. Molenda et al., Solid State Ionics 12, 473 (1989). 

[4] H. Yakabe et al., in Proc. 16th Int. Conf. Thermo- 
electrics, (IEEE, Piscataway, 1998) pp. 523-527. 

[5] I. Terasaki, Y. Sasago, and K. Uchinokura, Phys. Rev. B 
56, R12685 (1997). 

[6] R. Funahashi et al, Jpn. J. Appl. Phys. 39, L1127 (2000). 

[7] A. C. Masset et al. Phys. Rev. B 62, 166 (2000). 

[8] S. Hebert et al, Phys. Rev. B 64, 172101 (2001). 

[9] A. Maignan et al, J. Appl. Phys. 92, 1964 (2002). 
[10] D. Pelloquin et al, J. Solid State Chem. 170, 374 (2003). 
[11] T. Yamamoto, K. Uchinokura, and I. Tsukada, Phys. 

Rev. B 65, 184434 (2002). 
[12] Y. Miyazaki et al, Jpn. J. Appl. Phys. 41, L849 (2002). 



[13] J. Sugiyama et al, Phys. Rev. B 66, 134413 (2002). 
[14] J. Sugiyama et al, Phys. Rev. B 68, 134423 (2003). 
[15] J. Sugiyama et al, Phys. Rev. B 67, 214420 (2003). 
[16] J. Sugiyama et al, (2003) unpublished. 
[17] H. R. Krishnamurthy et al, Phys. Rev. Lett. 64, 950 
(1990). 

[18] M. Fujita, M. Ichimura and K. Nakao, J. Phys. Soc. 

Jpn. 60, 2831 (1991); M. Fujita, T. Nakanishi, and K. 

Machida, Phys. Rev. B 45, 2190 (1992). 
[19] M. Mikami, thesis, Osaka University (2003) in Japanese. 
[20] T. Tani et al, J. Mater. Chem. 13, 1865 (2003). 
[21] Y. J. Uemura, Muon Science edited by S. L. Lee et al, 

(Institute of Physics Publishing, Bristol, 1999) pp. 85- 

114. 

[22] G. M. Kalvius, D. R. Noakes, and O. Hartmann, Hand- 



■5 



book on the Physics and Chemistry of Rare Earths edited 
by K. A. Gschneidner Jr., L. Eyring, and G. H. Lander, 
(North-Holland, Amsterdam, 2001) vol. 32, chap. 206. 

[23] G. Griiner, Density Waves in Solids (Addison- Wesley- 
Longmans, Reading, 1994) chap. 4. 

[24] I. Tomeno, and M. Oguchi, J. Phys. Soc. Jpn. 67, 318 



(1998). 

[25] D. J. Singh, Phys. Rev. B 61, 13397 (2000). 
[26] W. Koshibae, K. Tsutsui, and S. Maekawa, Phys. Rev. B 
62, 6869 (2000). 



